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Abstract 

This paper is concerned with the well-posedness analysis of the Hartree-Fock 
system modeling the time evolution of a quantum system comprised of 
fermions. We consider quantum states with finite mass and finite kinetic 
energy, and the self-consistent potential is the unbounded Coulomb interac- 
tion. This model is first formulated as a semi-linear evolution problem for 
the one-particle density matrix operator lying in the space of Hermitian trace 
class operators. Using semigroup techniques and generalized Lieb-Thierring 
inequalities we then prove global existence and uniqueness of mild and clas- 
sical solutions. To this end we prove that the quadratic Hartree-Fock terms 
are locally Lipschitz in the space of trace class operators with finite kinetic 
energy. 

Technically, the main challenge stems from considering the model as an evo- 
lution problem for operators. Hence, many standard tools of PDE-analysis 
(density results, e.g.) are not readily available for the density matrix formal- 
ism. 
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I Introduction 



The time dependent Hartree-Fock theory provides approximate evolution 
models for many-body quantum systems comprised of fermions, as it ac- 
counts for the Pauli exclusion principle. It was first derived by Dirac pQ and 
simplified by Slater [2J. The Hartree-Fock system is a non-linear evolution 
equation for the one-particle density matrix operator g(t). It has the form 
of a von Neumann equation: 

lQ t {t) = [H(t),Q(t)], t>0 
0(0) = ? 

for a given initial value g 1 , and [., .] denotes the commutator of operators. The 
Hamiltonian H(t) = H(g(t)) is implicitly time- dependent, since the system is 
self-consistently coupled to the Poisson equation in this one-particle picture. 

A mixed state of a quantum system is usually described by a positive, Her- 
mitian trace class operator g(t), named density matrix operator, acting on 
L 2 (M 3 ) (see jlj). Hence, g(t) is usually represented as an integral operator 
with kernel g(x,y,t): 

(g(t)f)(x) = f g(x,y,t) f(y)dy V/a 2 (l 3 ). 

R3 

The Hamiltonian of the Hartree-Fock model can be written as 

H(t) = H + V H (t) - V HF (t), (1.2) 

with the operators 

H = — g A free Hamiltonian (without restriction of generality 

we set the Planck constant equal to one), 

V H = V H [g] Hartree potential, 

V HF = V HF [g\ Hartree-Fock potential. 

V H (t) is a (local) multiplication operator by the real valued function 

V H (x,t) = — [ ^\dz . (1.3) 

47T J \X - Z\ 

Here, n(z, t) denotes the (real valued) particle density of the system g(t) and, 
formally, it is obtained by n(z, t) = g(z, z, t). Note that the Hartree potential 
(II. 3J) is the Newtonian potential solution of the Poisson equation 

AV H = -n. (1.4) 



2 



The (non-local) Hartree-Fock-correction term (or exchange part) 

[V HF (t), g] = V HF o q — q o V HF is an integral operator whose first term 
has the kernel 

(kernel o §)) (x,y,t) = ± J <?(*, y, t) dz , (1.5) 

and hence V" HF [f3] has the kernel 

V HF {x,z,t) = — g^lfl^ . (1.6) 

47T |X — 2 1 

Here, we used 'o' to emphasize the composition of operators. 

Often it is convenient to rewrite the initial value problem (jl.ljl as an evolution 
problem (integro-differential equation) for the kernel g of g: 

Qt = H x g - H y g , t > 0, 

g(x,y,t=0) = g'(x,y), 

with g 1 denoting the kernel of g 1 . The subscripts x and y indicate that the 
Hamiltionian H acts, respectively, only on the x or the y variable: 

H x (t) = -±A x + V H (x,t)-V x HF (t) ■ (1.8) 

Here, V" F (t) is the integral operator (|1.5|) with kernel V HF (x, z, t). Analo- 
gously, V" F (t) has kernel V HF (z, y, t). The terms V H og and V HF og are obvi- 
ously quadratic in g which is the main challenge for an existence-uniqueness 
analysis of or (|1.7|) . On a first glance it would look easier to analyze 
the time evolution of the density matrix function g(x, y, t) according to (|1.7|) 
rather than the evolution of the operator g{t) . However, the main problem is 
to control the "diagonal" n(x,t) of g(x,y,t) without including (unphysically) 
many spatial derivations into the function space for g. As we shall see in Sec. 
HH this "control" of n(x, t) occurs very naturally when g(t) is a trace class 
operator (cf. [H], [0] for a more detailed discussion). The resulting draw-back 
is that the analysis of the operator evolution equation (jl.ljl is technically 
much more involved than analyzing an integro-differential equation of type 

O). 

We remark that there is generally a third approach for analyzing a Hamil- 
tonian quantum system of form (jl.ljl . A self-adjoint trace class operator g 
has a complete orthonormal system {yjj}j g pj C L 2 (IR 3 ) of eigenfunctions with 
corresponding eigenvalues {Aj} je N £ ^(N). Its operator kernel then has the 
"diagonal" representation 

Qfay) = Yl X i¥j(x)¥j(y) ■ (1-9) 
i&N 
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One easily verifies that, due to the Hamiltonian form of the eigenvalues 
\j of g(t) are constant in time (cf. |7|). Hence, can be rewritten as the 
following Schrodinger system for the evolution of the eigenfunctions ifj(x,t): 

i&<Pj(x,t) = H x (t)tpj(x,t), j'GN, t>0 
Vj (x,t = 0) = (p^x), JEN. 

Here, cp'j, Xj are the eigenfunctions and eigenvalues of g 1 . From ()1.9|) the 
particle density becomes 

n(x,t) = g(x,x,t) = ^2 ^il^iOM)! 2 , ( Ln ) 

and hence the Hartree-Fock Hamiltonian reads: 
H x (t)<pj(x, t) = --A<pj(x, t) 



+ J_ x f \<Pk(z,t)\ 2 <Pj(x,t) - <p k (x,t)<p k (z,t)(pj(z,t) dz 
An k J \x — z\ 

where the are a-priorly known from the diagonalization of g 1 . In (ll.K)jl 
the (fj evolve according to the same Hamiltonian H x (t), and they are only 
coupled through the Hartree and Hartree-Fock potential terms. 

The formulation (Jl.lUj) of the Hartree-Fock model is well suited for a mathe- 
matical analysis, since (jl.llj) yields a rigorous definition of the particle den- 
sity in L\R 3 ). Well-posedness of fpTTUJ) in H^R 3 ) and H 2 (R 3 ) was proved 
for Coulomb interactions in [S], and extended to more general interaction 
potentials in jH]. ^D] analyzes the corresponding Hartree model in L 2 (R 3 ), 
and the semiclassical limit and large-time behavior of (jl.lOj) is investigated 
in [H]. Further, quasiperiodic solutions to the Hartree-Fock system in an 
external electro- magnetic field are constructed in [T2] . 

We point out that is almost equivalent to (jl.lj) : the unique solution 

of also solves (jl.lj) . but uniqueness of the solution g{t) of (jl.lj) does 

not follow. Hence we shall directly analyze (jl.lj) in this paper, following 
the approach of Ref. 0. This strategy was also used in ^3] to prove exis- 
tence and uniqueness of a trace class operator-solution to the Hartree-Fock 
system. In that paper, however, the Coulombian two-particle interaction po- 
tential (appearing in (jl.3j) and (jl.8j) ) was approximated by a bounded 
function which greatly simplified the analysis (The Hartree-Fock term is then 
locally Lipschitz in the space of trace class operators, and the kinetic energy is 
not needed). A second motivation for our approach is that the reformulation 
(jl.lOj) becomes impossible for open quantum systems, since the eigenvalues 
Xj would then be time-dependent. Open quantum systems are important in 
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many fields of applications (quantum diffusion, coupling to a heat bath, cf. 
|14j . [To] . and references therein) and they are modelled by augmenting 
the right hand side of (jl.lj) by interaction terms of Lindblad form (cf. [T7j . 
[TS|). In our subsequent analysis we shall not include such (bounded) Lind- 
blad operators, but they would not pose any additional analytical problems. 

The paper is organized as follows: in Sec. |H]we introduce the functional set- 
ting for our subsequent analysis. Using perturbation techniques from semi- 
group theory we prove the existence of local-in-time solutions for in Sec. 
11111 In Sec. HVI we give a rigorous proof that our solutions to the Hartree-Fock 
system are mass and energy conserving. These a-priori estimates then imply 
that the constructed solutions are global. 



II Notations and functional setting 

We shall use the notation X x and X 2 for the spaces of, resp., trace class op- 
erators and Hilbert-Schmidt operators acting on L 2 (M 3 ). They are equipped 
with the norms (see also [H?] ) 

Mh ■= Tr\g\ , (2.1) 
Uh ■= (Tr|£| 2 ) 1/2 . (2.2) 

Since H = —A/2 > , the operator H^ 2 := ^/A is well-defined (having 
the symbol |£|/\/2 in Fourier space) and we can introduce the normed linear 
spaces 



Z := {g G X 1 | g is Hermitian, H^ 2 qHIJ 2 G 1±} 



and 



I nil ■- HUH _1_ III U-l/2-rrl/2 HI 

\Q\\Z ■— \\Q\\l + Ill-no 



Y := {g G I 1 \ g is Hermitian, H g G Xi}, 

\\q\\y ■= III ^ in i + \\H e\h 



(2.3) 



(2.4) 



Here, Hl^ 2 gH^ 2 denotes the closure (on L 2 (M 3 )) of the operator H^ 2 gH^ 2 , 
which is only defined on (a subset of) H 1 (R 3 ). Due to the compactness of 
q G Z, however, it can be extended to a bounded operator on L 2 (M 3 ). In 
the sequel we shall mostly suppress this closure symbols to keep the notation 
simple. In Lemma 12.11 we shall show that Y <Z Z holds . 

We shall denote operators in the form g (with an overwritten "hat") to dis- 
tinguish them from their kernels. For a self-adjoint operator g G X 2 its 
kernel g has the diagonal representation (jl.9j) . where {Aj}j g N G l 2 (N) and 
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the complete orthonormal system {fj}jeN C L 2 (M 3 ) are, resp., the eigenval- 
ues and eigenfunctions of g. g is Hermitian if and only if its kernel satisfies 
g(x,y) = g(y,x) (the bar denotes complex conjugation). For self-adjoint op- 
erators g G 1\ we even have {A-,}.,-^ £ Recalling the definition of the 
particle density in (jl.llj) . it is now possible to estimate n(x) via the trace 
norm of g: 

\\n\\ L i(R3) < l\g\\i = \Xj\ , 

with an equality for g positive. This natural control of the L 1 -norm of n 
constitutes the main justification for considering the Hartree-Fock system 
(HUD in X 1 instead of the PDE ljT77|) . 

While physical quantum states only lie in the cone of positive operators of Z 
or Y , we shall consider here the whole spaces as this simplifies the subsequent 
analysis. From a physical point of view the space Z comprises quantum states 
with finite mass (i.e. Trf? < oo) and finite kinetic energy, which is defined as 

EUQ) ■= ^{H]' 2 g H]' 2 ) ■ (2.5) 

in A 1/2 a A 1/2 

For g > this equals \\H gH Moreover, we can compute the kinetic 
energy (J2.5j) in terms of the eigenvector decomposition of g as 

EUlel) = \zZ HV^(x)||i 2(R3) = \\Hl' 2 eH l J 2 \\x (2.6) 

jeN 

(cf. Ref. 0, Lemma A.l). Hence, g G Z implies 

{^■} i6N C H\R 3 ), (2.7) 

and as a consequence 

.. 1 1 2 || 1 1 2 ii 1 1 2 



III Olll IT 1 / 2 ~ TJ 1 / 2 H 



l • 



(2.8) 

which is equivalent to the Z-norm. 

As we shall see in the next section, initial conditions in the spaces Z and Y 
give rise to, resp., mild and classical solutions of 

Some properties of Z and Y are stated in 
Lemma 2.1. 

(a) Z is a real Banach space. 

(b) Y is a real Banach space. 

(c) Y is a dense subspace of Z . 
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Proof, (a) The linear subspace of Hermitian trace class operators is closed 
in X x . Now let {Qj}jeN be a Cauchy sequence in Z. Hence 3 a, 7 G I 1 such 
that, for j — ► 00, 

& — <r, F 1/2 ^ /2 ^ 7 inZ x , (2.9) 

and the corresponding Cauchy sequences of the kernels satisfy, as j ; — ► 00: 

^■(a;, y) — > a(x, y) , H$H$ ej (x, y) — > y) in L 2 (M 6 ) . 

By Fourier transforming in both x and 7/ (and denoting its dual variables by 
£ and 77) we have 

— ^(£,77) inL 2 (M| xMj,(l + -^\ v \ 2 )didr,) , 
which is a complete space. Therefore 

l(x,y) = H^H^ y 2 a(x,y) and gj — ► a in Z. 

(b) Same argument as for (a) (see Ref. fHj for details). 

(c) The proof is divided into two steps: first we define an auxiliary space 
X cY, then we prove that X is dense in Z. Let 

X := {g E Ji I g Hermitian, AgA E XJ, (2.10) 

\\g\\ x := Uh + ll A £ A llli ■ (2-H) 

iStep 1: show X CY. 

Let pel and decompose it (and at the same time AgA) in its positive and 
negative parts: Q = Q+ — Q-, AgA = Ag + A—Ag_A with g±>0,Ag±A> 
0. From g±,Ag±A E I\ follows y/]}±, Aa/^± G X 2 . Hence A£± G Xi, which 
implies g E Y . 

Step 2: show that X is dense in Z. 

Let g E Z . Assume without restriction of generality that g > (otherwise 
separate into g±). We have 

V 

n^iu = y^^j(iiyjiii 2 (R3) + 2ii v ^iii 2 (K 3 ))' 

where Aj, ^ are the eigenvalues and eigenfunctions of g. For all e > : 
3iV G N with ^ satisfying \\g — g N \\z < e/2, with the kernel g N (x,y) = 
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Y^j=i^jfj( x )fj(v)- For eac h n G N we consider approximations a n , which 
we define in terms of their kernels: 

N 
3=1 

with appropriate functions ip™ G H 2 (K 3 ) such that ip™ — ™ ipj in H 1 (M 3 ). 
Note, that we do not require {ip^}^ =l to be orthonormal. We have a n G X 
(since G H 2 (R 3 ) and rank(er n ) < N), o n > as sum of positive operators. 
Further, 

a n ™ g N and H^a n H^ ™ H^g N H^ 
in the strong operator topology, since we have V/ G L 2 (IR 3 ): 

N 



(*»/)(*) = E W*) I m^(y)dy 

N r 

™ / /(i/)vi(y)dy = (ff Ar /)(*)mi 2 (» 8 ), 



N 



(H^a n H^f)(x) ^XjH^ix) [ f(y)H^(y)dy 



{Hl' 2 Q N H l J 2 f)(x) in L 2 (R 3 ). 



Their norms satisfy 

N N 

III ^"n III i = Tr<x n = Aj||^" |||2( R3 ) — > E^'ll^'lli 2 (K 3 ) = 1^ I 1 ' 

iV 

_^ > A ; I I \ { I ' ~ 

A 



i4 1/2 ^o 1/2 ii = ^E A iii v ^ii^ 



|E^ii v ^-ii^») = i^o /2 ^o 1/a Bi- 



2 

i=i 

The two equalities on the left hand side are easily verified by evaluating the 
trace in an arbitrary orthonormal system of L 2 (IR 3 ). By Grumm's conver- 
gence theorem ([20 , Theorem 2.19) we obtain 

III " -AT in ™-™ n III 171/2,2. £Vl/2 rrl/2 ~N rrl/2 1|| n~>°° n 

iFn - £ 111 — ► U, \\\H ' a n n > — tl Q ' g n Q ' ||i — > U . 
Since a n rt —^ g N ' m z, X is dense in Z. And so is Y. □ 
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Ill Local-In-Time Solution 



This section is concerned with the local-in-time solution of the Hartree-Fock 
system We rewrite as 

gt(t) = -i[H ,g{t)} - i[V H (t),g(t)} + i[V HF (t),g(t)} 

= h (g) + F(g) , (3.1) 
g(0) = Q 1 , 

where we formally define the operators h and F as: 



h (g) := -i [H , g(t)) , 

F(g) := -i \V H {t) - V HF (t),g(t)} ■ 

First we consider the free evolution equation 

Qt(t) = h (g) , t > 0, 

Q(0) = Q 1 



(3.2) 



(3.3) 



in the "energy spaces" Z and Y . This linear problem admits a unique global 
solution in X\ (cf. [22] Chap. 5). It can be represented via the isometric 
Co— evolution group {G (t),t e R}, which reads 

G (t)g = e~ i6ot g e i6ot . (3.4) 
Its infinitesimal generator h is defined as 

T>(ho) = {g e 2i | gV(H ) C V(H ) = if 2 (R 3 ), (H g - gH ) is an operator 
with domain H 2 (R 3 ) with an L 2 (R 3 ) extension H g — gH G Xi} , 
h (g) = -i[H ,g]. 

We remark that Y C V(h ), while Z £ V(h ). The group {G a (t), t E R} 
preserves Hermit icity and posit ivity in Xx(cf. Ref. [22] Theorem 5.1). 

The following Theorem 13. II states that the restriction of {G (t),t G R} to Z 
and Y yields a global solution of ()3.3|) in these two spaces. 



Theorem 3.1. The evolution group G restricted to Z (resp. Y) is an iso- 
metric C — evolution group on Z (resp. Y). 

* 1/2 

Proof. Since H and G (t) commute, it follows directly from the corre- 
sponding properties of G (t) on X x . □ 
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As a second step we shall consider ([3.1)1 as a perturbation of ([3.3)1 by the 
perturbation operator F (see j2I] Theorem 6.1.4). We shall demonstrate that 
F is locally Lipschitz in Z and in Y, which guarantees the existence of a 
unique local- in-time solution for 1)3.1)1 . as stated in Theorem 13.61 
Next we show that the perturbation F is locally Lipschitz. It is demonstrated 
by Corollary 13.51 which is the result of the following three lemmata where 
we prove that F maps Z into Z, and resp., Y into Y. 

In the sequel C denotes generic, but not necessarily equal constants. will 
denote the weak LP— space (cf. [23J, e.g.). 

Lemma 3.2. Let geZ. Then V HF [g] G B(L 2 (R 3 )). 

Proof. The eigenf unctions of g satisfy (cf. p. 9)1 . 1)2.7)1 ) <p 3 - G H 1 (R 3 ), which is 
continuously embedded in L P (M 3 ) for 2 < p < 6. For arbitrary / G L 2 (M 3 ) we 
use (fjf G L p (R 3 ) for 1 < p < 3/2 and l/\x\ G L 3 (M 3 ). Then the generalized 
Young inequality yields (where V denotes the convolution operator): 



\x\ 



! ) , 3 < q < oo , 



(3.5) 



and the Holder inequality implies: 



G L P (I 



6 

- < p < 6 
5 



(3.6) 



Using the kernel p. 6)1 of V" HF we estimate for all / G L 2 (IR 3 ): 

e(x,z) 



^\\V HF [g]f\\L^) 



\x — 2 



/(*) dz 



L 2 



< 



3& 



< Cbllz||/|| L2 



\x — z\ 



f(z) dz 



L 2 



| 


W 


L 3 (R 3 ) 




x)f(x)) * 


1 

\x\ 






L 3 (R 3 ) 




L 6 / 5 (R 3 ) 


1 



L 6 (R 3 ) 



£2. 



X 



i2 



l/l 



L 2 



where we used f)2.8)l for the last estimate. Hence, ||V ahf [^]||b(l 2 (r 3 )) < C||f?IU- 



□ 



10 



Lemma 3.3. (a) Let g G Z. Then Hi 12 V hf o g Hi' 2 G 1 x . 
(b) Let geY. Then H V HF ogel,. 

" 1/2 

Proof. For the subsequent analysis we first note that H and V are "equiv- 
alent operators" in the sense that 

3 

Hi' 2 = K i d *i and 9 Xj = -V2HWRJ. (3.7) 

3=1 

Here Kj and Rj (j = 1,2,3) are bounded pseudo-differential operators de- 
fined by their symbols (in Fourier space) 

Rj denotes the j-th Riesz transform operator in M. 3 (cf. |24| . p. 58 for defi- 
nition and details). Since kj,rj G L°°(]R 3 ) we have Kj,Rj G B(L 2 (R 3 )), and 
hence: 

Hi' 2 g Hi 12 el, if and only if d Xj g d Xk G 1, \/j,k = 1,2,3. 



Part (a): To prove the assertion we will show 

d Xk V HF o g d x . el, Vj,k = 1,2,3. 

Since V HF is an integral operator with kernel ()1.6|) . 9 Xfc y HF is also an integral 
operator with the kernel 

x T /H Ff \ 1 d x k Q(x,y) 1 g(x,y) 

d Xk V (x,y) = j : — | (x k - y k ) . 

Air \x — y\ Air \x — y\ 6 

The last two terms define, resp., the functions D*(x,y) and —D*(x,y). For 

the two integral operators D\ and belonging to these kernels we will now 
show: 

D*gd Xj el, Vj, k =1,2,3, (3.9) 

D k 2 gd Xj e 1, Vj, k =1,2,3, (3.10) 
which will then imply 

d Xk V HF ogd Xj = b\gd Xj + D\ g d Xj G 1, V j, k =1,2,3. 

To prove ()3.9|) we consider the identity 

d* e a,, = (#y 2 + 1)" 1 (#y 2 e + g d Xj ) . (3.11) 
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In (13.11)1 . HlJ 2 Q d x . G I] follows with (|3.7[) immediately from g G Z. To 
prove g d Xj G 1 1 we decompose the self-adjoint operator g into its positive 
and its negative parts: g = g + — g_ with g± > 0. Hence we have 



Q± Hp = g¥ (e W) e X, 



since both factors of the right hand side are in X 2 (for the second factor we 
have Hl*Q±Hl* G X x if and only if q^h\* G X 2 , cf. Ref. jH], Lemma A.l). 
(13.7)1 then yields £ c^. G X x and it remains to show that t)\ (Hq + I) G 
B(L 2 (R 3 )). Since (i/ 1/2 + /)~ 1 is a bounded operator from L 2 (R 3 ) to iX^R 3 ) 

we set for any / G L 2 (M 3 ) : [h]' 2 + f)" 1 / =: g G fl^R 3 ). Hence it suffices 
to prove 

|| A* #IU 2 (» 3 ) < C Ibll^CR') G H\R 3 ) . 
Using the eigenfunctions representation ()1.9|) of g's integral kernel we have 



4nD*g{z) = X j n - rM) fP^Kg(y)dy. (3.12) 

Since (pj G H 1 (R 3 ) for g £ Z (see (12.7)) ) it remains to prove that the second 
(integral) factor of (|3.12|) is in L°°(IR 3 ). Since we proceed here like in the 
proof of Lemma l3~^l we only give the key estimates. By a Sobolev embedding 
we have 

<p d , g G H\R 3 ) ^ L P (R 3 ) , 2<p<6. 
Hence, ()3.5|) directly gives 

(jfj ' 9) * 1 — r e L q (R 3 ) , 3 < q < 00 . 
Next we consider its spatial derivatives: 

(d Xk Tp- y g + Jp-. d Xk a)*-^ e I/(M 3 ) , 3 < p < 00 

by proceeding as before. In total we have [Tpj ■ g) * l/\x\ G H /1,P (1R 3 ) , 3 < 
p < 00 , and the following estimate holds due to a Sobolev embedding: 

N'^*rL( K 3) - C NLncR-) 11^11^3) ■ ( 3 - 13 ) 
Now we can estimate ()3.12|) : 

|| D* p|| L2(R 3) < C l^'l H^illi 1 ^ 3 ) Hfi'll^MR 3 ) 

< c \\q\\z \\g\\m 
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This proves D k {Hi 12 + I) 1 G £(L 2 (M 3 )), and hence £ <9 2 . G X x . 

The proof of (|3.1U|) is analogous to the proof of (J3.9)) : In order to show 
D k (H}J 2 + I)" 1 G £(L 2 (R 3 )), we consider 



4nD k 2 g(x) = J2 \ I ' pr^Ts ( x h ~ Vk) 9(v) dy, V # G if 1 (M 3 ). (3. 14) 

i6N / 3 |X 

To estimate the last integral we use x k /\x\ z G L 3 J 2 (R 3 ) andTpJ-g G L P (R 3 ) , 1 < 
p < 3, and get by the generalized Young inequality 



x k 



\x\ 



3 



Finally, 



47r ||^2 ^|| L 2(R3) ^ ^ l A jl ll^ilUW (<Pj-9)* 

< C7 HffHz |MI^(B") 



L6(R3) 



which completes the proof of Part (a). 

Part (b): By the same technique as before we shall prove d 2 x V HF o g G 
Xi, k = 1,2, 3. We consider the kernel of the integral operator d 2 V HF : 



(3.16) 

where the right hand side defines the kernels of three integral operators A\, 
A\, A k (in this order). To prove Afg G I±, k,l = 1,2,3, we shall use the 
factorization 

A k g = A?(H + I)- l (H g + g), 

with H g + g G T\ since g G Y. Hence, we have to show that Af (H + J) -1 G 
B(L 2 (R 3 )), that is 

\\A k g\\ L ^) < C\\g\\ H2m Wf e L 2 (R 3 ), 
where g G H 2 (R 3 ) is defined as # := (H + J) -1 / . 

Case 1: Boundedness of A k (H + I) -1 . In order to deal with the second 
derivatives of g(x,y) in A^, we introduce two self-adjoint operators: 

Si := -{H g + gH ) G X 1; 6 2 := ^(H g - gH ) G X 1; 
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with their kernels satisfying 61,82 G L 2 (IR 6 ). The eigenvector decomposition 
of Si and 5i yields 

^i{x,y) = ^2fij^j(x)%() j (y) and 5 2 (x, y) = <TjXj(x)xj(y), (3-17) 

jeN jeN 



with 



IIWIIi 



\\H g + gH \\i < l\H gl\i 



|# £ - ^i7 ||l < |||# £ 



II > 



(3.18) 
(3.19) 



and {ipj}, {Xj} are orthonormal systems in L 2 (IR 3 ). For each g G if 2 (R 3 ) 
L q (M?), 2 < q < 00, we consider 



fj ■= (ifjjg) * — e L P (M 3 ), 3 < p < 00 and 

3> 



A Sobolev embedding then implies 

||/.j|U°°(R3) < C||V'il|L 2 (R3)||5'||H2(IR3) . 

The same argument holds for hj := (x]g) * l/\ x I : 

IIMl°°(r 3 ) < c|lxilU 2 (R 3 ) Ibllff^R 3 ) • 

Using H g = 5i + i5 2 and (|3.17|) - (|3.21j) we estimate: 



(3.20) 



(3.21) 



A x g(x,y) 



g{y)dy 



jeN 



\x - y\ 



\x - y\ 
dy + i^2a jXj (x) J 



jeN 



Xj(y)g(y) 

\x - y\ 



dy 



L 2 (R3) 



< Yl 1^1 n^iiU 2 (R3)ii/iiu-(R3) + n iixiiU 2 (R 3 )ii^-iu 

jeN jeN 

< cJ2(\^\ + WM\g\\HH^) < 

jeN 



00 1 



tf 2 (R3) < C||£||r||fl'||ff2(K3) • 



Case 2: Boundedness of A^Hq + 1) 1 . With the eigenfunction decompo- 
sition of g we have 

A\g = -7T Y] Aj d Xk <fj(x) / ^-(y) n: ^ 0(3/) °^ 

27r ~^ JR3 
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As in (|3.15|) we have (ipjg) * (x k /\x\ 3 ) G L P (M 3 ), 3/2 < p < oo, and its spatial 
derivatives are in L q (M. 3 ), 3/2 < q < 6. Hence we have 



< C||^||hi(r3)||^| 



So we obtain the desired result 

II^2#IU 2 (r 3 ) < c l^'l II^jIIhvk 3 )!!^!!^ 2 ^ 3 ) - cil^lklls'll^ 2 

Case 5: Boundedness of A^(H + 1)" 1 . We rewrite the last term of (|3.16j) as 



jeN 1 1 



and estimate as in Case 2: 



< C1lf?l|y|Mltf 2 



and hence 

||^3fl'||L 2 (]R3) < C ^ |Aj| ||yj||l,6(lR3) 

- IlljilnW) ll^l|g 2 ( 

Summarising we have proved H V HF (H +I)~ 1 G i3(L 2 (IR 3 )) and the assertion 
(b) follows. □ 

Lemma 3.4. F(g) is Hermitian for a given g G Z . 

Proof. Let g G Z, hence it is Hermitian. The assertion is a simple conse- 
quence of 



V H [g](x) G R, V HF [g](x,z) = V HF [g](z,x) 
i.e. the self-adjointness of V H [g] and V HF [g] (cf. (JS3)). 



□ 



The analogous properties of the previous Lemmata for V^[£>] can be proved 
with the following generalization of the Lieb-Thirring inequalities (cf. Refs. 
|S] and Theorem A. 3.): 

IMIl, ( r3) < C Mil {MHl /2 \g\Hl /2 )) l - a , i<g<3, 

||Vn|U. (K3) < C r \\gfx {^{Hl'XeWl 12 )) 1 '^ 1 < r < 3/2, 
a = (3 - q)/2q, = (3 - 2r)/2r, 

As a consequence of the Lemmata 13.21 - 13.41 we then obtain: 
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Corollary 3.5. F maps Z (resp. Y) into itself and is locally Lipschitz in Z 
(resp. Y). 

Using standard perturbation results (Ref. [21], Theorem 6.4.1,) the local 
Lipschitz continuity of the function F(g) guarantees the existence of a local- 
in-time solution of (|3.ip . 

Theorem 3.6. Let g 1 G Z. 

(a) Then the Hartree-Fock system \3. 1)) has a unique mild solution g G 
C([0,t max );Z) with a potential V H - V HF G C([0, t max ); B(L 2 (R 3 ))). More- 
over, if the maximum time interval is finite, i.e. t max < oo, then 

limtTW. llfflU = 00 ■ 

(b) In case g 1 G Y C V{h Q ), g is a classical solution with g G C([0, t max ); Y)fl 
C\[0,t 

mux 

(c) For a// < ti < t max i/ie map ^ 7 i— > g(t) is Lipschitz continuous on some 
(small enough) ball {\\g — g \\z < ^(^i)} C Z , uniformly for < t < t\. 

Remark 3.7. Since equation (|3.1|> is in Hamiltonian form, it preserves pos- 
itivity: if the initial value g 1 > 0, then g(t) > V£ G [0,t maa; ). 

From now on we will restrict our analysis to positive density matrix op- 
erators, which represent physical quantum states. 

IV Global-In-Time Solution 

In order to prove the global-in-time existence of solutions to equation (J3.ll) . 
we shall derive an a priori estimate for the kinetic energy (cf. Lemma f4.2|) . 
This estimate is a consequence of the conservative character of the prob- 
lem. More precisely, we show that the total charge and the total energy are 
conserved by the local-in-time solution g of the Hartree-Fock system. On 
a formal level this is well known since Dirac [1J, but we shall need here a 
rigorous proof. We recall that the total charge corresponds, by definition, to 
the quantity \\g\li. And the total energy of the Hartree-Fock system is given 

by 

Etot{o) '■= Ekin(g) + E pot (g) , 

where the kinetic energy Eki n (g) is defined in (|2.5j) and the potential energy 
equals 

E po t(Q) ■= \t±(qV*) - ~Tr(#n . 
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Lemma 4.1. Let g 1 G Z and g 1 > in hS. Then the local solution g of 
Theorem VJ.tA satisfies 

\W)\\i = Wh, E tot (g(t)) = Et^g 1 ) Vt G [0,t max ) . 

Proof. The conservation of the total charge for (J3.1|) can be proved by writing 
the integral equation associated to (|3.1|) in terms of the semigroup {G (s)} 
and the perturbation F(g): 

Q(t) =G (t)g I + [ G (t-s)F(g(s))ds, < t < t max (4.1) 
Jo 

and taking the trace. For g(s) G Z we have V H (s) - V HF (s) G B(L 2 (R 3 )). 
Using the representation (|3.4|) of G (t) and the properties of the trace we 
obtain TrF(g(s)) = 0, which implies Trg(t) = Trg 1 Vt G [0,t max ). 

Concerning the conservation of the total energy we first restrict the analysis 
to g 1 G Y, which by Theorem 13.61 implies the local-in-time existence of a 
classical solution g(t). The general case g 1 G Z will be then derived by a 
density argument. For g G Y all of the following manipulations are well 
defined. Hence we apply the operator H to both sides of ()4.1|) and then take 
traces. Commuting G (t) with H , using V H (s) - V HF (s) G B(L 2 (R 3 )), and 
using the cyclicity of the trace, we obtain (as an extension of Lemma 3.7 in 
Ref. 0): 

Tr(H g(t)) = Tr(#,#) - % [ Tr(H V H [g} g(s) - H () gV H [g](s)) ds 
t Jo 

+i [ Tr(H V HF [g] g(s) - H {) gV HF [g}(s)) ds 

t 
= TriH^+if Tr(V H {g] [H , g{s)\) ds 

t 
-i [ Tr(V HF [g] [H ,g(s)}) ds 

Jo f t 

= Tr(F ^)- / Tr(V H ig}g t (s)) ds+ [ Ti(V HF [g] g t (s)) ds, 
Jo Jo 

(4.2) 

where we used ()3.1|) in the last step. Applying an integration by parts in t 
(possible since g t G C([0, t max \X-^), the second term on the right hand side 
of ()4.2|) becomes 

Tr(U H [£] g t (s)) ds 

-Tr(V H [g](0) g 1 ) + Tr(V H [g](t) g(t)) ~ f Tr(V t H [g](s) ds. 

Jo 

(4.3) 
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For this last integral term we note that n t G C([0, t max ); /^(R 3 )) implies, 
by the Hardy-Littlewood-Sobolev inequality, that V t H [g] = (1/4tt|z|) * n t G 
C([0,t max ); L 3 (IR 3 )). Using g G C([0,t m a X );Y) and techniques like in the 
proof of Lemma l3~3l one easily obtains V^f^p G C([0, t max );Ii) and the 
following equivalence holds by using the Poisson equation (jl.4|) and two in- 
tegration by parts (in x): 

Tr(V t H [g] g(s)) = - [ V t H [g](x, s) AV H [g](x, s) dx 

Jm. 3 



1 d 
2dt 
1 d 
2di 
1 d 
2d~t 



\VV H [g}{x,s)\ 2 dx 



V h [q](x,s) n(x,s) dx 
Tr(V s [g]g(s)) . (4.4) 



Since g t G C([0,t ma x);Tx) and V HF [g] G B{L 2 {R 3 )) (cf. Lemma the last 
term of (|4.2|) is well defined. Using ()1.5|) we finally compute it: 

Tr(V HF [g] ft(a)) = -!-(/ dx / dz f^'f ft(*,s,a)) 



47r 2dtVy R3 J R 3 [x — z| 

= ~Tr(V"*\e] fa)) . (4.5) 

Plugging the expressions (|4.3|) - (|4.5jl into (|4.2jl we obtain 

Tr(tf ^)) + ~Tr(t^[£] g(t)) ~ \t±(V*'[q) g(t)) 

= Trfag 1 ) + ~TY(V^](0) g 1 ) - ^Tr(U^](0) &) , 
which proves the assertion for g 1 G F. 



Next we consider an initial condition g 1 G Z with g 1 > 0. Due to Lemma 
I2.1f c) it can be approximated by a sequence {g' n } C F such that g' n — ► g 1 , 
with ^ > 0, Vn G N. For any < ti < t max , the corresponding trajectories 
Qn(t), < t < t\ lie in F (cf. Theorem I3.6f b)) and they converge in Z to 
£?(t), uniformly on [0, t\] (cf. Theorem I3.6f c)). This implies convergence of 
the kinetic energy of g n = g n (t) Wt G [0, ti\: 

E km (Qn) = \\H l J 2 g n Hl' 2 h ™ l^e^lli = E Ue) ■ 
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To prove convergence of the Hartree-Fock potential energy we estimate: 
\Tr(V HF [g n ]og n )-Tr(V HF [g]og)\ 

< \\V HF [g n } ~ V HF [g}\\ B{L 2 (R3)) \lg~n\lx + Hy^telllB^CR 8 ))!^ ~ eh 

< C\\g\\z \\Qn-Q\\z n ^ 0, (4.6) 

where we used the estimate ||V AHF [^ n — q] \\b(l 2 (r 3 )) < C\\g n — q\\z obtained in 
the proof of Lemma 13.21 

The convergence of the Hartree potential energy, i.e. 

is obtained analogously by using the estimate ||^ H [£] ||l°°(m 3 ) < C||£IU de- 
rived in Equation (3.46) of Ref. [5]. 
The assertion of the lemma now follows from 

E kin {g n {t)) + E pot {g n {t)) = E km {g' n ) + £ poi ($J, Vn 6 N, t E [0, h] 

in the limit n —>■ oo. □ 

Lemma 4.2. Let g 1 e Z and g 1 > 0. Then the kinetic energy of the Hartree- 
Fock system \3. 1\) is bounded V< G [0, t max ): 

< E kin (g(t)) < E tot {g') . (4.7) 

Proof. The kinetic energy Eki n {g{t)) is non negative by definition and equals 
Etotig 1 ) — E pot (g(t)) by Lemma l4~2l To complete the proof we have to show 
that E pot (g(t)) > 0, Vt E \0,t max ). 

2E pot (g(t)) = Tr(gV H )-Tr(gV HF ) 

= £Aj / dx\^( X )\U]- [ dz ^ 



x — z\ 



47r|x — z| J K 3 J]R3 47r|x — z 

which is non-negative since, by the Schwarz inequality, we have: 
\q(%,z)\ 2 = I > A:'" tyAx) X]'~ ibAz 



* ( 

= n(x) n(z) . 
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□ 



From Lemma \4 . 1 1 and l4~2*l we conclude that the local- in-time solution of The- 
orem can be extended to t max = oo and the global solution theorem 
follows: 

Theorem 4.3. Let g 1 G Z and g 1 > 0. Then the Hartree-Fock system 
\3. 1)) has a unique global mild solution g G C([0, oo); Z) with a potential 
V H — V HF g C([0, oo); i3(L 2 (M 3 ))). 

In case g 1 G F and g 1 > 0, g is a global classical solution with g G C([0, oo); Y)f] 
r'( 0.^):I). 
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